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Abstract 

Given a Heath-Jarrow-Morton (HJM) interest rate model M and a 
parametrized family of finite dimensional forward rate curves Q, this pa- 
per provides a technique for projecting the infinite dimensional forward 
rate curve rt given by M onto the finite dimensional manifold QJYhe 
Stratonovich dynamics of the projected finite dimensional forward curve 
are derived and it is shown that, under the regularity conditions, the given 
Stratonovich differential equation has a unique strong solution. Moreover, 
this projection leads to an efficient algorithm for implicit parametric es- 
timation of the infinite dimensional HJM model. The feasibility of this 
method is demonstrated by applying the generalized method of moments. 



1 Introduction 

Consider a complete filtered probability space (0, T, (-Ft)j P) satisfying the usual 
conditions, and where P denotes the risk-neutral measure. Consider an HJM 
model M. (see [5]) with the Musiela parametrization (see JJ)- The dynamics of 
the forward rate curve under P can be given by the following infinite dimensional 
stochastic differential equation (SDE): 

dr t (x) = p,(r t ,x)dt + a(r t ,x)dW t , (1.1) 

where W is an m-dimensional standard P-Brownian motion. From the arbitrage- 
free condition given by [5], it follows that 

d f x 
fl(r t ,x) = — r t (x) +a(r t ,x) J a(r t ,u) T du. 
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Since is an infinite dimensional stochastic differential equation, it is 

difficult in practice to apply it to characterize the dynamics of the forward rate. 
More commonly, simpler parametric families of curves are used to characterize 
the forward rate curve. This practice gives rise to the following issues: 

• Under what conditions, if any, does the infinite dimensional forward rate 
curve r t (-) defined in Hl.lfl admit a finite-dimensional realization? 

• Under what conditions is this forward curve consistent with a given family 
of finite-dimensional curves? 

Under the assumption of smoothness of the volatility function a, Bjork et al. 
(12,0) give necessary and sufficient conditions for the consistency problem, and 
for the existence of a finite-dimensional realization. Their results are developed 
further by Filipovic and Teichmann ([7]). (See also Filipovic JBJ and Bayraktar 
et. al. 1 .) In this paper however, we will assume that the forward rate curve 
does not admit a finite dimensional realization. It is of interest to find the pro- 
jection of the infinite dimensional forward rate curve onto a finite dimensional 
manifold, and we consider this problem here. Precisely, starting from the speci- 
fied finite dimensional manifold, each time the velocity vector points out of the 
this manifold, it is projected onto the tangent space. The velocity vector field 
obtained by this procedure now describes a flow which is finite dimensional. 

When the finite dimensional manifold is linear this projection characterizes 
the behavior of constructing forward rate curves from cross sectional market 
data by interpolation using a parameterized family of finite dimensional curves. 
Using our results this data can be used to calibrate an infinite dimensional 
HJM model (see also Remark l3.2[) by deriving the coordinate dynamics of the 
projected forward curve provides an efficient way for implicitly estimating the 
evolution of the original HJM model. 

The remainder of the paper is organized as follows. In Section 2, we in- 
troduce a fundamental mathematical formalism within which we consider this 
problem. In particular, we confine our discussion to the case in which the solu- 
tion r t (-) given by lies in a specific Sobolev space with a finite time horizon. 
Then in Section 3, we derive the dynamics of the coordinates of the projected 
forward curve given a finite dimensional manifold Q. In Section 4, we propose 
a simple and efficient strategy for parametric estimation of the original HJM 
model through the projected forward curves. In particular, we demonstrate 
that the generalized method of moments can be straightforwardly implemented 
in this situation. 
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2 Basic Mathematical Formalism and Projec- 
tion 



To ensure sufficient regularity of model (|l.lfl . we restrict attention to the situ- 
ation in which the solution r t (-) lie within a Sobolev space B with a finite time 
horizon To defined as the space of all infinitely differentiable functions 

r : [0, T ] ^ R 

having finite norm 

°° r T ° ( d n r \ 2 

Remark 2.1 In order to include all constant curves into the space, there usu- 
ally exists an exponential weighting junction e~ lx ( r y > 0) in the definition of 
the norm for a weighted Sobolev space. However, in our case, since we already 
confine us in a finite time horizon, this weighting function becomes unnecessary. 

Now let us consider the forward rate dynamics r t (•) satisfying Equation ■ 
By using the infinite dimensional Ito formula (see can be written in 

Stratonovich form as 

dr t (x)= (^{r t ,x)-^K(r t ){a{r t ))]{x)^dt + a{r t ,x)odW u (2.1) 
where o denotes the Stratonovich integral and 

m 

o-' r (r t )a{r t ) :=^a' ir {r t )(Ti{rt). 

i=l 

Here a' r denotes the Frechet derivative of the vector field a with respect to r t . 
For convenience, let us denote 

M (r t ) := A(n) - ^.(r t >(r t ). (2-2) 

The following assumption guarantees the existence a local strong solution for 
Hl.lj l in the space B. 

Assumption 2.1 i) a is a smooth mapping; 

ii) the mapping r i— * er(r) J Q o~(r,u) T du — ^a' r (r)(a(r)) is also smooth. 

By a smooth mapping, we mean this mapping belongs to C°°. 
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Remark 2.2 From Assumvtion \2.1\ it is straightforward to deduce that 



fi(r,-)£B, a{r,-)eB, VreB. 

Therefore, within a regular space B, we can start to discuss the projection of 
the forward rate flow r t onto a certain finite dimensional manifold. 

Let Q :— {G(z, •) € B : z = (z 1 , z n ) 6 Z C K"} be a finite dimensional 
manifold, where Z is the parameter space. Consider a curve j3 in this manifold of 
the form (3 : h i— > G{z{h) 1 •), where ft. i— > is an n-dimensional curve. Then, 
by the chain rule, we can compute the following Frechet derivative: 



Dp(0) = DG(z(h)r)\h=o 



k=l 



dz 



Therefore, for each zq in Z, the tangent space in this manifold is then ob- 
tained by considering all the curves passing through G(zq, •) and tangential to 
Q. From this argument it is clear that the tangent vector space is given by 



f dG(z,-) 



dGjz, ■ 



Without loss of generality, it is assumed that dG S Z i '^ . 



are linearly in- 



dependent. 

Let us now consider orthogonal projections between any linear space V € B 
and the tangent vector space S% . Since the generating basis of S? is not 
necessarily orthogonal, we have the following projection formula II^ : V <— > 

9z" / ' 



f 9G(z,-) 

span | -±r^ 



E 

i=l 



dGjz, 
dzi 



dG(z, 



(2.3) 



where (Ay) := A" 1 with A := ((^ d ,^T i )j- Here (•,•} denotes the inner 
product of the projection. It is worth mentioning that this inner product is 
not necessarily the one defined in the space B, which was constructed so as to 
guarantee the boundedness of d/dx, together with Assumption H2.1fKiD assures 
the smoothness of fi. Actually an appropriate choice of this inner product should 
reflect the calibration criteria of a financial institution. Generally we can define 
a Hilbert space TL with a finite time horizon Tq containing all the functions 



v:[0,T o ] 
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satisfying \\v\\n < oo with norm || • \ \-h denned as 



\Mn ■= / v 2 (x)w(x)dx, (2.4) 
Jo 

where w(-) is a positive bounded continuous function on [0, To]. The correspond- 
ing inner product is defined by 

rTo 

(u,v) = / u(x)v(x)w(x)dx, \/u,vEH. (2-5) 
Jo 

Typically, fund managers, when interpolating the forward curve by using a 
certain curve family would require more accuracy in the long term region, and 
accordingly would choose choose w(x) to give higher weight to larger values of 
x. Alternatively, traders would choose w(x) to put greater weight on the short 
term. 

Remark 2.3 It follows from that B C H. 

Assumption 2.2 We assume that G(z, •) G B. 

Remark 2.4 By Remark \2.!& and Assumvtion XH.'lA we have that 

a(G(z t ,-))eB, n(G(2t),-)€B 
dG(z t ,-) 

and — — - — G n, V % = 1, 2, n. 

oz 1 

Therefore, we can discuss the projection without worrying about its existence. 
The projected forward rate flow of rt onto a finite manifold Q is thus given 

by 

dG(z u •) = IT- \ji(G(z t , •))] dt + n 2 - t [a(G(z u •))] ° dW t . (2.6) 
where II^ is defined by (|2.3(l . 

3 The Dynamics of Finite Dimensional Coordi- 
nates 

In order to derive the diffusion process for z, first we give the following lemma. 

Lemma 3.1 For any two continuous semimartingales X t and Y t having finite 
quadratic variation on finite intervals, i.e. [X]t < oo, and \Y]t < oo for each 
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T > 0. Then for V T > 0, we have 



(X t Y t )odW t = X t o(Y t odW t ), a.8. 



(3.1) 



Proof. From the definitions of Stratonovich's integral and variation, it is suffi- 
cient to prove that 

71 

lim y y (X ti -X u _ 1 )(Y ti -Y u _ 1 ){W u -W u _ 1 ) = 0, a.s., (3.2) 

where "J = {to,ti, ...,t n } is a partition of [0, T] and ||^|| is the mesh of this 
partition, namely, ||A|| = maxi<i<„{|ii — Since for any A, we have 



Y^{x u -x ti _j(y ti -y^m-w^) 



< max \W U - W u _,\x 
1< i <r 



\ «=1 



(3.3) 



by the uniform continuity of the Brownian motion W on the compact support 
[0,T] and the fmiteness of [X]t, and [Y]t, it follows that i|3.2|) is true. Lemma 
follows. □ 
By applying the Stratonovich chain rule, we can derive that 



dz l 



(3.4) 



On the other hand, by l|2.3[) and l|2.6[) . we also have 



dG(z t , 



E 

i=l 



Emmg^,-)) 



E 

i=l 



^A y (a(G(z t ,-)), 
i=i 



dG(z t: 



dGjzu ■) 
dz* 

dG(z t , ■ 



cH 



dz 



odW t (3.5) 
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By comparing (|3.4|) and t\'d. 5fl and using Lemma I3~T1 we conclude that the finite 
dimensional vector z follows a diffusion process: 



dz\ = 



3=1 



dG{z u - 



^A y -(a(G(z t ,-)) 

3 = 1 



dt 



odW t , V i e {!,..., n}.(3.6) 



Theorem 3.1 Suppose G satisfies Assumvtion \2.'2\ and that Assumvtion \2. 1\ is 
valid. Then there exists a unique strong solution for the SDE fff.6J) . 

Proof. This result follows directly from Remark |2. 21 and Proposition 5.2.21 in 

ma. 

Remark 3.1 i?ac/i trading day banks perform a yield curve fitting to determine 
the prices of bonds with maturities whose prices are not observed because they 
are not traded on the market. This is done because the payoffs of even vanilla 
products such as caps, floors and swaptions consists of a sequence of cash flows. 
And the valuation of several products like these require the entire yield curve. 
Each trading day the prices of several bonds are set by the yield curve obtained by 
curve fitting. Assume now that the prices between the two calibrations follow an 
HJM model (this models the market activity between the times of curve fitting). 
We will show the calibration of the banks can be obtained by discretizing \2.b]) . 
i.e., we will show that G{zAti •) given by \2.b)) in fact minimizes 

min \\r A t ~ G(yAt,x)\\-H- (3.7) 
G(j?A*,-)ee 

Therefore using the results of next section one can calibrate an infinite dimen- 
sional HJM model using the yield curve data obtained form weigted linear re- 
gression. 

The minimization in \3. 7| ) is equivalent to showing 

(r(At) - G(z A r),G z (z A t, -)) H = 0, (3.8) 

where G z is the derivative of G with respect to z. Assume that r S Q and take 
let the SDE in \2. b}) also start from this point. Then discretizing \2.1)) we obtain 

rAt = ro + n(r ,x)At + a(r ,x)WAt- 



7 



We can apply the same discretization to \2. b)) and obtain 

r A t - G(zau •) = Ofa, x) - n 2 - Q \p(G(z , •))]) At + (a(x) - n 2 - n [a(x)]) W At . 

(3.9) 

By the definition of the projection operator ILg the right-hand- side is perpen- 
dicular to the S? therefore so is the left-hand- side, hence $3.8p is satisfied. 

Remark 3.2 Consider the following affine term structure: 

G(z t ,x)=g (x)+ gi (x)zl + ...+g n (x)z? (3.10) 

Then the minimizer in jff.TI ) is guaranteed to be unique, since Q is linear. 

4 Parametric Estimation of the Original HJM 
Models 

In this section, we discuss parametric estimation of the original HJM model M. 
by using its projection onto a finite dimension manifold Q. Generally speaking, 
estimating parameters of an infinite dimensional SDE is a difficult task to im- 
plement. However, since we have already derived the SDE i|3.6[l governing the 
dynamics of its projected finite dimensional vector z, this provides us an easier 
way to empirically investigate the underlying HJM model. 

Suppose 9 is a vector of parameters parametrizing the HJM model. By 
projecting this forward curve into a finite dimensional manifold Q, we have the 
following general form for the diffusion process of z: 

dz t =A(z u 8) + B(z t ,9)odW t , (4.1) 

where W is an m-dimensional standard Brownian motion. It will be convenient 
to write l|4.1|l in Ito form which is as follows. 

dz t = (a(z, 9) + l - J2 (?§K^)) B j (z t ,9)\ dt + B{z u 9)dW u (4.2) 

where B^ denotes the jth column of B. 

Here we apply the generalized method of moments (GMM) proposed by 
Hansen (1982 [5]) to estimate 9. For banks, periodic calibration of the initial 
forward curve produces a time series {z tk }i<k<N- Let A := tk+\ — tk, for each 
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k = 1, N. By discretizing 1)4. we obtain a discrete-time model: 



(z tk ,0))B^z tk ,0) A + B(z tk ,0)e k , 



(4.3) 

for k — 1,2, ...,N, where e k is a m-dimensional Gaussian random vector with 
mean and covariance matrix JA, i.e., {efc}i<fc<Ar are mutually independent. 
Since z tk and e k are independent we construct the moment functions {h k (0)} 

as 



h k {9) = z tk+1 - z tk - I A(z tk ,9) + 2 2. ( — ) Bi(z tk ,6) | A, (4.4) 



for k — 1,...,N and denote the sample average by /at(#) :— -A> Xl&Li hk{9)- 
Then by 

N = min{(f N (0)J N (0))}, (4.5) 

8 

is the least squares estimator of 8. Under fairly general conditions (see [H]), the 
estimator On offers a consistent estimator of 0q. 

Remark 4.1 // the dimension of is high, it is straightforward to strengthen 
this algorithm by adding more moment functions. 

As argued by Hansen in the estimator l|4.5|l is generally not efficient as far 
as its convergence rate is concerned. The least squares estimator of 1(4.5(1 can be 
improved by taking a weighted least squares estimator. Suppose that {h k (0o)} 
is strictly stationary, where 9q is the true parameter value, and define 

T„ := E {h k (9 )h k - v (9 y} . (4.6) 
Assuming {r„} is absolutely summable we define 

oo 

S:= ]T IV (4.7) 

U — — OO 

The optimal GMM estimator is given by 

0* N = min{f N (0)S- 1 f^(0)}. (4.8) 

u 

To implement (14.8(1 we need an inital estimate of S, which calls for an initial 
estimate 0. The least-squares estimate 1(4. 5|l gives such an estimate, which can 



9 



be used to calculate an estimate for S as follows (flU): 

q 

v,N + T^AT ) (4.9) 

I/=l 

where 

1 N 

Fv ' N= N Yl h n {6)h n - v {0)' (4.10) 

and where 9 is the estimator obtained using (|4.5|l . This estimate of S is then 
used to compute an estimate of 9 using Q4.8[) . This obvious recursion can then 
be carried until the estimator becomes stable. 
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